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i> ■ ABSTRACT 

a^ 

The discovery of Gamma-Ray Burst (GRB) afterglows revealed a new class of 
Q I variable sources at optical and radio wavelengths. At present, the environment and 

Q I precise redshift of the detected afterglows are still unknown. We show that if a GRB 

^f) ' source resides in a compact (^ 100 pc) gas-rich environment, the afterglow spectrum 

will show time-dependent absorption features due to the gradual ionization of the 
surrounding medium by the afterglow radiation. Detection of this time-dependence 
can be used to constrain the size and density of the surrounding gaseous system. For 
f~^ , example, the Mgll absorption line detected in GRB970508 should have weakened 

^N I considerably during the first month if the absorption occurred in a gas cloud of size 

j~^ I ^ 100 pc around the source. The time-dependent HI or metal absorption features 

^^ ' provide a precise determination of the GRB redshift. 

^ ■ 
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1. Introduction 

The detection of delayed emission in X-ray (Costa et al. 1997), optical (van Paradijs et 
al. 1997, Bond 1997), and radio (Frail et al. 1997) wavelengths following Gamma-Ray Bursts 
(GRBs), the so-called "afterglows", can be reasonably explained by models in which the bursts 
are produced by relativistically expanding fireballs (Paczyhski & Rhoads 1993; Meszaros & 
Rees 1997; Vietri 1997a; Waxman 1997a,b; Wijers, Rees, & Meszaros 1997; Vietri 1997b; Sari 
1997). On encountering an external medium, the relativistic shell which emitted the initial 
GRB decelerates and converts its bulk kinetic energy to synchrotron radiation, giving rise to the 
afterglow. Since the emission occurs on time and length scales which are much larger than those 
of the initial explosion, the dynamics of the expanding shell is self-similar and thus insensitive 
to the properties of the energy source (Blandford & McKee 1976). The source size inferred 
from radio scintillation data of GRB970508 (Frail et al. 1997; Taylor et al. 1997) is consistent 
with the generic prediction of this model (Goodman 1997; Waxman, Kulkarni, & Frail 1997). 
Moreover, synchrotron emission by Fermi-accelerated electrons behind the expanding shock can 
successfully explain the gross properties of the time-dependent spectra observed in the afterglows 



of GRB970228 and GRB970508 for a fireball energy ~ 10^^ gj.g (Waxman 1997a,b; Wijers et al. 
1997; Sari, Piran, & Narayan 1997). 

In the simplest emission model, the afterglow flux is proportional to the square root of the 
ambient medium density, and could therefore span a dynamic range of more than five decades if 
the ambient density varies between its characteristic value in dense molecular clouds (~ 10^ cm~^) 
and in the intergalactic medium (~ 10~^ cm~^). Even merely galactic environments span a 
wide range of possible gas densities, between those of ellipticals and spirals. The possibility that 
GRB sources reside in both gas-rich and gas-poor environments might account for the absence of 
afterglow emission in some GRBs (e.g., Groot et al. 1997). Any direct observational verification 
of the postulated connection between the afterglow flux and the GRB environment would provide 
a crucial test of the basic fireball model. 

Despite the expectation that bright afterglow sources should reside in gas-rich environments, 
there is currently no direct measurement of the density, size, or redshift of the host system for 
known afterglows. The discovery of absorption lines in the optical afterglow of GRB970508 
provided the first direct estimate of source distance, constraining its redshift to the range 
0.835 < z <^ 2.3 (Metzger, et al. 1997; Djorgovski et al. 1997). The possible existence of host 
galaxies has yet to be conclusively determined from direct imaging of afterglow fields (Sahu et al. 
1997; Fruchter et al. 1997). 

If a GRB source resides within the HI disk of a spiral galaxy, its optical-UV spectrum is 
likely to show a Lyman-limit trough, a damped Lya absorption, and metal absorption lines, all 
at the source redshift. However, the same absorption features could also be produced by an 
intervening galaxy along the line-of-sight to the source. The ambiguity between the GRB host and 
an interloper could be removed by detecting the Lya forest which should end at the true source 
redshift. However, for source redshift -z ^ 2, such a spectroscopic detection cannot be made with 
ground-based telescopes and requires the use of the Hubble Space Telescope (hence being feasible 
only for bright afterglows, due to the small collecting area). In this Letter we argue that if the 
absorption occurs in a sufficiently compact region (^ 100 pc) around the GRB source itself, then 
the afterglow radiation will ionize the region and change the observed absorption features with 
time. The release of ~ 3 x 10^^ ergs in ionizing photons could ionize ~ 10^^ hydrogen atoms (or 
~ IO^Mq), and create an ionized bubble of radius ~ 100 pc n^ in the surrounding galaxy, where 
ni is the ambient proton density in cm~^. Given the typical HI column densities found in galactic 
disks (Broeils & van Woerden 1994), such a bubble could be sufficiently big so as to poke a hole in 
the disk and eliminate the absorption features of the host altogether. Aside from establishing the 
GRB redshift, detection of this time-dependent effect could be used to constrain the gas density 
and size of the GRB host. GRB afterglows are much more effective at ionizing their environments 
than supernovae, because their peak UV luminosity is higher by several orders of magnitude. 

The computation of the time-dependence of the absorption spectrum is a non-trivial 



problem^ A snapshot of the ionizing effect that the radiation shell has on the ambient medium 
is schematically illustrated in Figure 1. At the late times under consideration, the blast wave 
that produced this radiation is already non-relativistic and lags behind the radiation front. Only 
the front of the afterglow light pulse moves into the unperturbed gas. Each of the subsequent 
radiation shells propagates into a medium which had already been partially ionized by prior 
shells. Since each radiation shell has a different spectrum at emission and a flux that declines 
with radius as l/r^, the ionization state of the ambient medium is generally inhomogeneous and 
time-dependent. The eventual absorption spectrum seen by a distant observer, can be calculated 
by coupling the radiative transfer equation to the radiative ionization equation for the medium. 
In §2, we formulate the set of causal equations that describe the processed afterglow spectrum, 
and in §3, we present numerical results for the time evolution of the associated absorption lines. 
Our discussion ignores recombination, since the recombination time (~ 10^ yr nj~ ), is typically 
much longer than the lifetime of the afterglow. Our main conclusions are summarized in §4. 



2. Ionization Dynamics 

We consider a GRB source which turns on at time t = and illuminates a stationary ambient 
medium of uniform density n, with a time dependent luminosity per unit frequency, Lu(t). Since 
most of the interstellar absorption occurs at large radii, we ignore the initial sub-parsec fireball 
phase during which the radiation is still coupled to the hydrodynamics. In the late phase of 
interest, the pressure wave has already slowed down to sub-relativistic speeds, so that it lags 
behind the radiation field and has no effect on the absorption spectrum of the outgoing radiation 
shell. The ionization state of the gas is then affected only by the luminosity Li^{t) at frequencies 
higher than the ionization threshold z^q- 

The photoionization rate per atom of species a (=hydrogen, helium, or metals) is given by, 

1 dnair,t) f^ F^{r,t) 

- du — (yu,a{y) , (1) 



na{r,t) dt Juo,a hu 

where Fy{r,t) = Ly{i)/ATTr'^ is the flux at a radius r from the source, and o"bf(i^) is the 
photoionization (bound-free) cross-section at a photon frequency u. All photons are assumed 
to propagate radially away from the central source. Equation (||) determines the evolution of 
the number density of atoms, na{t), from that of the flux, Fu{r,t). This flux is absorbed in the 
medium through both bound-free and bound-bound transitions. The radiative transfer equation 
in spherical symmetry reads, 

l_ d[r'Fu{r,t)] ^ _^^(^^^)^^^(^^^) [cTM.aM +^bb,aH] , (2) 



j,2 ^j. 



Note that the observed time, which is related to the width of the outgoing radiation shell, it typically much 
shorter than the light crossing time of the absorber. 



where Cbb is the cross-section for bound-bound absorption. 

The coupled equations (|l|) and (||) determine the evolution and spatial distribution of the 

flux and the ionization fraction, for a given set of initial and boundary conditions. We assume for 

simplicity that the source is embedded in a spherical gas cloud of uniform density and radius R, 

namely 

t. n \ i "-".O for r < i? 
Uait = 0,r) = < , (3) 

I otherwise 

and start the calculation at a radius tq ~ 10^ cm with the boundary condition, 

FAn,.t) = !^. (4) 

Our results are insensitive to the precise choice of tq, as long as this radius is much smaller than 
the scale over which the ambient interstellar medium is distributed. 

Following Waxman ( 1997a, b), we model the time and frequency dependence of the afterglow 
luminosity at emission as 



where, 

'^mit) = 1.7 X 10- (i|)' {^y^EliXr Hz . (6) 

Here ^b and ^e are the fractions of the equipartition energy in magnetic field and accelerated 
electrons, E = 10^'^E^2 erg is the fireball energy, iju- = (i/hr), and 

L.^ = 8.65 X 10^9 V^ (^] E,2 -^, (7) 



0.1/ sec Hz' 

where ni is the ambient proton density in units of 1 cm^^. In our numerical examples we use the 
typical parameter values which are required to fit the existing afterglow data, namely ^£-£'52 = 0-2, 
^B = 0.1, -y/ni-E'52 = 1, a = 0.5 for u > Um, and a = —1/3 for v < Vm- 

Finally, the atomic cross-sections in equation (g) are calculated as follows. 

(i) Bound-Bound absorption. The cross-section for transitions between two bound states, i and k, 
of an atom has a Lorentzian shape (Rybicky & Lightman 1979) 

where Xk,i = c/i'k,i is the transition wavelength, and A^^j is the transition probability which is 
related to the absorption oscillator strength /^ j through 

Ak,^ = 6.670 X 10^5 gg^-i 5fc/^,t/[5fc(Afc,,/A)2] , 



with gi and gk being the statistical weights of the levels i and k. The values of these parameters 
for permitted resonance lines of the most important astrophysical elements are summarized by 
Verner, Verner & Ferland (1996). To minimize the number of free parameters, we assume a cold 
quiescent absorber. Thermal or turbulent broadening would change the line profile but not affect 
its time dependence, which is dictated by bound-free transitions [cf. Eq. (|l])]. 

(ii) Bound-Free absorption. An analytical expression for the photoionization cross section exists 
for single-electron atoms (Hall 1936), and reads for hydrogen, 

/ X A /i^o\'^exp[4- (4tan-ie)/e] „ ^ ,^, 

\ ly J 1 — exp (— ZTT/e) 

where Aq = 6.30 x 10~^^ cm^ and e = ^/u/uq - 1, with uq = 3.29 x 10^^ Hz being the 
photoionization frequency threshold for an electron in the ground state. For simplicity, we assume 
that all neutral hydrogen atoms are in the ground state; the excitation rate to higher levels does 
not exceed the photoionization rate for the characteristic power-law spectrum of GRB afterglows. 

For non-hydrogenic atoms, Verner et al. (1993) obtained a useful analytical fit to the 
photoionization cross section ct„^/(z^) for an electron in a shell with quantum numbers n,l. The 
cross section can be written as an^i{i^) = a^, F{y) for frequencies v > Vn,h where y = v jv^. Here 
Vy^^i is the threshold frequency for photoionization of an electron in the shell n, /; a* = 0">,(n, /, Z, N') 
and Vi, = i'^,{n,l, Z,N) are fitting parameters which depend on the atomic number Z and number 

/ s -P 

of electrons N; and F{y) = [{y — 1)^ + y^] y~^ (l -|- ^y/ya) , where Q = 5.5 + I — 0.5P, and 
Vw, Va, and P are three additional fit parameters which depend on n, I, Z and A^. The total cross 
section to photoionization is obtained by summing over all shells, 

'7bf(l^) = ^(^n,i{i^) ■ (10) 

n,l 

We ignore secondary ionization by fast electrons which are created by X-ray or 7-ray 
photoionization. 



3. Numerical Results: Time Dependence of Absorption Lines 

Since hydrogen is by far the most abundant element in galaxies, we examine first its prominent 
Lya absorption signature. Figure 2 shows the expected evolution of the afterglow spectrum in a 
narrow frequency band around the Lya resonance for an HI column density of 3 x lO^'^ cm~^ and 
two different cases of spherical absorbing clouds around the source: (i) a compact {R = 1 pc) and 
dense (nH,o = 100 cm"'^) cloud (left panels), and (ii) a more extended {R = 100 pc) and rarefied 
(™H,o = 1 cm~^) cloud (right panels). The first case illustrates the conditions in a molecular cloud, 
while the second corresponds to a galactic disk of scale height ~ 100 pc. Although the two spectra 
possess the same absorption profile initially due to their identical column density, their time 



evolution is very different. The smaller absorber is ionized faster, because its atoms are exposed to 
a higher ionizing flux and because it requires fewer photons to get ionized due to its smaller mass. 

The photoionization efficiency can be quantified through the change in the equivalent width 
(EW) of the absorption line, defined as 

EW= [dX (^'-■0-^'-) , (11) 

where Fu{X) is the observed spectral flux (at r ^ i?) as a function of wavelength across the line 
and F^^q{X) is the continuum flux that the afterglow would have had if there were no absorption. 
The values of EW for the Lya line are listed in the different panels of Figure 2. For a compact 
absorbing region such as a molecular cloud, the EW could decrease by orders of magnitude within 
several days after the GRB detection. 

Although ground based observations are not capable of detecting Lya absorption for absorber 
redshifts z ^ 2, they can still probe metal absorption lines as demonstrated in the case of 
GRB970508 (Metzger et al. 1997; Djorgovski et al. 1997). Because metal lines are often 
unsaturated, it is easier to detect the time evolution in their case than for the damped Lya trough. 
Figure 3 illustrates the expected evolution of the Mgll absorption line at Aq = 2798 A, which was 
detected in the spectrum of GRB970508. Our calculations assume a solar abundance of Mgll 
(Anders & Grevesse 1989). Since there are far fewer Mgll atoms than HI atoms, it takes fewer 
photons and less time from the onset of the afterglow, to ionize them. Therefore, the evolution of 
the metal line EW in Figure 3 is more rapid than that of the Lya EW in Figure 2. 

In Figure 4 we show a broader range of the spectrum, including the UV and X-ray regimes, 
where flux is absorbed by photoionization of hydrogen. For a ~ Ipc absorber, the absorption 
profile narrows considerably after a couple of weeks. Consequently, the observed flux from a source 
at Zs = I in the emission band of 13.6-150eV increases from 10~^^ erg cm~^ s~^ after 6 hr to 
4.2 X 10~^^ erg cm~^ s~^ after 2 days, despite the overall fading of the afterglow emission. 

All of the above time dependences occur only for absorbers which are associated with the 
GRB source. More distant absorbers along the line-of-sight are exposed to an afterglow flux which 
is much too weak to affect their ionization state during the brief period of time over which the 
radiation lasts. 

If the GRB environment is metal rich, then it is likely to contain also dust. The broad-band 
extinction by dust would suppress the observed optical- UV flux from the afterglow by a wavelength 
dependent factor. The deviation from the fireball predictions would weaken at late times due to the 
heating of dust above its sublimation temperature by the afterglow radiation (Waxman &: Phinney 
1997); this process might explain the transient suppression of optical relative to X-ray emission 
during the first two days of GRB970508. The absence of an optical afterglow for GRB970828 
might have also been caused by dust extinction (Groot et al. 1997). This interpretation is 
supported by the low-energy turn-over in the ASCA X-ray spectrum of GRB970828, which implies 
a high hydrogen column density, ~ iq21-22 qt^-'^^ along the line-of-sight to the source (Murakami 



et al 1997). Figure 4 indicates that if this absorption originated from a compact molecular cloud 
around the GRB source, then the X-ray absorption feature must have evolved on a timescale of 
days. For a given inferred column density, the level of evolution becomes weaker with increasing 
absorber mass. 



4. Conclusions 

Afterglow models postulate the existence of gas-rich environments around some GRB sources. 
We have shown that this conjecture can be tested empirically by monitoring the afterglow 
absorption spectra. If the surrounding gas is not fully-ionized initially, then the bright afterglow 
radiation will modify its ionization state. The medium into which the afterglow radiation 
propagates at late times would be more ionized than the medium into which the early radiation 
propagated. Hence, the equivalent width (EW) of the associated absorption lines would decline 
with time during the afterglow (see Fig. 2). The rate of EW variation can be used to constrain 
the size of the absorbing region; the more compact the absorbing region around the GRB source 
is, the higher is the photon flux to which it is exposed, and the faster it gets ionized. When 
combined with the value of the column density inferred from the EW measurement itself, the 
size estimate can be used to find the density of the ambient medium^. One could then test 
statistically the predicted correlation between the afterglow flux and this density [cf. Eq. (^)]. In 
principle, it might even be possible to infer the profile of a non-uniform density distribution along 
the line-of-sight based on frequent monitoring of the equivalent-width evolution. Identification 
of time-dependence in the afterglow spectrum would isolate the variable (and hence, associated) 
absorption lines from other lines which are caused by chance intersections of unrelated absorbers 
along the line-of-sight, and hence provide a precise determination of the GRB redshift. 

Our results imply that if the observed metal absorption lines of GRB970508 (Metzger et al. 
1997; Djorgovski et al. 1997) occurred in a region of radius ^ 100 pc around the GRB source 
itself, then the equivalent width of the lines would have declined noticeably over the first month of 
the afterglow (Fig. 3). Similarly, the soft X-ray absorption which was possibly detected by ASCA 
in GRB970828 (Murakami et al. 1997) should have decayed under the same conditions (Fig. 4). 
Detection of these effects in future GRBs would help determine the properties of their gaseous 
environment and their redshift. 

Since the peak optical fluxes of afterglows such as GRB970508 are ^ 10^ times brighter than 
supernovae, they can be used to probe galactic environments at exceedingly high redshifts. If 
the GRB source population traces the star formation history of the Universe, then the dimmest 
GRB events known might originate from a redshift as high as z ~ 6 (Wijers et al. 1997). In 



^Note that if the absorber size or density can be measured independently by another method, then it might be 
possible to infer the luminosity distance to the source and by that constrain the values of cosmological parameters. 



hierarchical cold dark matter models for structure formation, high redshift hosts are expected to 
be denser [n oc (1 + zj)^] and more compact [R oc M^'^{1 + Zf)~^, where Zf is their (higher) 
formation redshift, and M is their (smaller) mass] than low redshift galaxies, hence making the 
effects discussed in this paper even more pronounced. 
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Fig. 1. — Snapshot of the outward propagating radiation shell and its effect on the ambient 
absorbing medium. The ambient gas is gradually ionized by the passing UV radiation. The column 
of absorbing atoms seen by the front of the radiation shell is higher than that seen by the later 
radiation. As a result, the equivalent width of the absorption lines declines with time from the 
point of view of a distant observer. 



11 



Fig. 2. — Evolution of the observed Lya absorption profile for a spherical gas cloud around a 
GRB source, with a radius R=l pc and an initial HI density nH,o = 100 cni^^ (left panels), or with 
i?=100 pc and nH,o = 1 cni^^ (right panels). The equivalent width (EW) of the line is listed in each 
case. The vertical axis shows the observed flux for a source redshift Zg = 1 (for an h = 0.5, fi = 1 
cosmology), and the horizontal axis refers to the emission wavelength at the source, A^ = A/(1 + Zs). 
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Fig. 3. — Evolution of the observed Mgll absorption line for a spherical gas cloud of solar metallicity 
with R=W pc and an initial hydrogen density of nH,o = 10 cm~^, or equivalently a Mgll density 
^Mg,o = 3.8 X 10~^cm~^ (left panels), as compared to a cloud with i?=100 pc, nu,o = 1 cm~^, and 
^Mg,o = 3.8 X 10~^cm~^ (right panels). Axis notations are the same as in Fig. 2. 
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Fig. 4. — Evolution of the soft X-ray and UV absorption for a spherical gas cloud with R=l pc 
and an initial HI density nn.o = 100 cm~^. The dashed line shows the afterglow flux that would 
have been observed if there were no absorption. Axis notations are the same as in Fig. 2. 



